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Abstract 



In this paper we treat the multiparticle hopping asymmetric diffusion model (MADM) 
on Z introduced by Sasamoto and Wadati in 1998. The transition probability of the 
MADM with N particles is provided by using the Bethe ansatz. The transition proba- 
bility is expressed as the sum of ./V-dimensional contour integrals of which contours are 
circles centered at the origin with restrictions on their radii. By using the transition 
probability we find P(x m (t) = x), the probability that the mth particle from the left is 
at x at time t. The probability ¥(x m (t) = x) is expressed as the sum of l^l -dimensional 
C*~) . contour integrals over all S C {1, • • • , N} with |5| > m, and is used to give the cur- 

rent distribution of the system. The mapping between the MADM and the pushing 
asymmetric simple exclusion process (PushASEP) is discussed. 



b ' 1 Introduction 

03 



The Bethe ansatz, which is useful to find the eigenvalues and eigenvectors of the Hamil- 
tonian of one-dimensional quantum spin systems, has been one of the main techniques in 
studying one-dimensional stochastic particle models [5ll6tl9HT6] . One direction of studying 
these models by the Bethe ansatz is the asymptotic analysis of the current that can be 
computed from transition probabilities. In this direction the simple exclusion process on Z 
has been extensively studied [6,8.ill U16l - [20] and some results are now used to show that the 
Kadar-Parisi-Zhang (KPZ) equation belongs to the KPZ universality class [3]. 

The model we are going to consider in this paper is the one- dimensional asymmetric 



diffusion model without exclusion. This model was also named the multiparticle hopping 
asymmetric diffusion model (MADM) by Sasamoto and Wadati [13]. The dynamics of the 
MADM with N particles is as follows: Each site on Z is equipped with 2N Poisson clocks 
with rates pr n and ql n (p + q = 1, n = 1, • • • , N) and all clocks are independent. If a clock 
with rate pr n {ql n ) rings at x, which is occupied by at least n particles, then n particles 
simultaneously jump to x + 1 (x — 1). However, if the number of particle at a site is less 
than n when the clock with rate pr n {ql n ) rings at the site, then nothing happens and the 
clock resumes. The rates are not free but are given by the (/-binomial coefficients 
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with g = j, where A + fj, = 1 and A, /u 7^ 0. These requirements on rates are needed for the 
Bethe ansatz solvability [H0Q3]. I n [H] only one free parameter was involved in the rates 
but it is possible to extend to have two free parameters p and A as above [2]. In [2] the 
model with two parameters has exclusion property and the multiparticle hopping in p3] 
is interpreted as pushing effect. In this pushing version the model becomes the drop-push 
model in the limit q, A — > [13] or more generally the pushing asymmetric simple exclusion 
process (PushASEP jl| in the limit A — > in [4]. 

The totally asymmetric simple exclusion process (TASEP) and the one-sided PushASEP 
are determinants in the sense that their transition probabilities are expressed as determi- 
nants whose entries are contour integrals [¥ [|15j . These transition probabilities are a starting 
point for the multi-point joint distributions as well as the currents of systems [Hl8|HH[T2]. 
In the asymmetric simple exclusion process (ASEP) with N particles, the transition proba- 
bility is not in a determinantal form but it is expressed as a sum of A^-dimensional contour 
integrals over all permutations in a symmetric group §jv [IB]- Hence, some techniques for 
determinantal models to study the current of the system are not applicable to the ASEP. 
In spite of this limitation, Tracy and Widom obtained the current distribution of the ASEP 
with step initial condition and found the fact that its fluctuation is governed by the GUE 
Tracy- Widom distribution [16H20] . 



The main goal of this paper is to take a first step to see the MADM is also governed 
by the Tracy- Widom distribution. As the first step, we provide the transition probability 
and the current distribution of the MADM with iV particles. In fact, the MADM has been 
expected to belong to the KPZ universality class by studying its energy gap [14] and it is 



1 In the PushASEP in A\, the pushing can occur in only one direction, so we will call it the one-sided 
PushASEP 



known that the current fluctuation of the one-sided PushASEP (equivalently the MADM 
in the limit A — > 0) with a special initial condition is governed by the Tracy- Widom distri- 
bution [2]. Hence, we expect our results in this paper to provide a starting point to extend 
a result of the one-sided PushASEP to the two-sided PushASEP (equivalent to the MADM). 

This paper is organized as follows. In Section [2] we find the transition probability of the 
MADM by using the Bethe ansatz. As in the ASEP |16| . the transition probability is given 
by the sum of A r -dimensional contour integrals. However, unlike the transition probability of 
the ASEP in [16] . there is a delicate condition on contours in the MADM. This difference is 
due to the difference in the S-matrices of the ASEP and the MADM, and the difference also 
affects the way of obtaining the current distribution. In Section [3] we provide the probability 
that the mth leftmost particle is at x at time t in Theorem 13.71 and the current distribution 
of the system in f|3.35j) . In Section H] we discuss the relation between the MADM and the 
two-sided PushASEP. The mapping between the MADM and the PushASEP is confirmed 
by their transition probabilities, and the probability that the mth leftmost particle in the 
two-sided PushASEP is at x at time t is provided. 



2 The transition probability of the MADM 



The MADM with N particles is a continuous-time Markov process with a countable state 
space. If we denote the mth leftmost particle's position by x m , then a state of the process 
can be specified by particles' positions X = (x±, ■ ■ ■ , x/v) € ^ N with x\ < ■ ■ ■ < xn- We 
will call the state space 

N 



{(xi,--- ,xn) eZ :xi<---<xat} 



the physical region of the MADM. By the standard technique |14H16| to use Bethe ansatz 
we have a differential equation for a function u(X; t) = u(x\, • • • , xjv; £) on Z x [0, oo) 



— u(X;t) = y^ \pu(xi, ■■■ , Xj-i, Xi - 1, x i+ i, ■■■ , x N ; t) (2.3) 

+ qu(xi,--- ,Xi-i,Xi + l,x i+ i,--- ,x N ;t) - u(X;t) 

with the condition imposed on the boundary of the physical region 

u(zi,--- ,X{,Xi -1, ••• ,XN',i) (2.4) 

= (m(xi,- ■■ ,Xi-l,Xi-l,--- ,x N ;t) + Xu(xi,--- ,Xi,Xi,--- ,x N ;t) 

for i = 1, • ■ ■ , N, and the initial condition in the physical region, that is, 

u(X; 0) = 8 Y (X) for x\ < ■ ■ ■ < x N , (2.5) 



where Y = (y 1} • • • ,y N ) is the initial state of the MADM. The solution of (|23|) with (|2T3|) 
and (|2.5p is the transition probability and we will denote it by Py (X; t). Define the S'-matrix 
for the MADM 



M + A^aC/3 



and let 



A„ 



n * 



ff(i)<r(j) 



(2.6) 



<r(i)>ff(j) 



for the coefficient of the Bethe ansatz solution of (|2,3p with (|2,4p . Here, <r is a permutation 
in a symmetric group Sat on {1, • • • , N} and the product is over all S-matrices that satisfies 
the given condition under the product symbol. If there is no S'-matrix that satisfies the 
condition, then A a is defined to be 1q Let 

em = ^ + qii~l, 

which is interpreted as energy of the MADM. 

Lemma 2.1. Let A + jj, = 1 (1/2 < [i < 1) and a < (3. Let C^ be a circle centered at 
with radius Ri in C and £j € C^ (i = a, /3). If 1 < R a < Rp < j, then 



Rp< 



§8-m 



Kp 



Proof. It is easy to see that 



Rp< 



Rp - \i 
XRb 



< 



£,p- V 



Kp 



where the first inequality is equivalent to 



1< Rr < V 



■p 



X 



D 



Lemma 2.2. Let A + jj, = 1 (1/2 < /i < 1) and a < (3. Let C^ be a circle centered at 
with radius Ri in C and £j € Cf^ (i = a, /3). If 1 < R a < Rp < c where 

I oo if X = 0, 



then 



Rr, < 



CaC/3 



A* + A^aC/3 



(2.7) 



2 In general, for any product of S-matrices with a certain condition, if there is no S'-matrix that satisfies 
the condition, then the product is defined to be f in this paper. 



Proof. Let A ^ 0. It is easy to see 



RgR/3 RgR/3 < 

/x + XRpRp /i + \R a Rp 
where the first inequality is equivalent to 

/' 



£a£/3 



fJ- + A^a^/3 



1< Rr < t-. 



£ 



A 



If A = 0, (HZD is obvious. 



D 



Lemma 2.3. Let X + fi = 1 (1/2 < /j, < 1) and a < (3 < 7. Let C^ 6e a circle centered at 
with radius Ri in C and & € C^ (z = a, /3, 7). 7/ 1 < 7? a < Rp < 7? 7 < c where 

U i/A^O 



then 



Rrv < 



I 00 if X = 0, 



/' 



(2i 



Proof. Let A 7^ 0. It is easy to see that 



R a Rj(l — XRj) R a R~f(l — XRp) 
K a < < < 



where the first inequality is equivalent to 



M 



£067 — -^oC/36) 



1< 7? 7 < ^. 

A 



If A = 0, d2SD is obvious. 
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Now, we find the transition probability of the MADM. 

Theorem 2.4. Let A + /x = 1 (1/2 < /i < 1) and Cj^ (i = 1, • ■ ■ ,N) be a circle oriented 
counterclockwise, centered at with radius Ri. Assume that 1 < R\ < • • • < Rn < c where 

c= U if A^O 

\oo if A = 0. 

The transition probability of the MADM with N > 2 particles is 



N 



w e 



dfi •••Cat. 



(2.9) 



Remark 2.1. In this paper we consider the nontrivial MADM with at least two particles 
because the MADM is just the simple random walk when N = 1. It can be easily shown 
that (|2.9p satisfies (|2.3[) and (j2.4|) for all X & 7j N in the same way as in the proof of the 
ASEP [16]. Hence, we shall prove only that (f2Uj) satisfies (f23j) . 



Proof. The initial condition (|2.5p is satisfied if 



which is the integral for the identity permutation id. Hence we will show that 

a^id JC RN JC Rl i 

for any state in the physical region. We prove (|2.10p by induction. First, we show that 

M + A66 - ft gi-K.-ig,-*-!^^ = (2 . n) 

when xi < X2 and yi < yi- Substituting ft = f- so that 77 runs over Cr where R = R1R2, 
the left hand side of (|2.1ip becomes 

/' + A '> - ^1 Tjn-ifi-iff-'i+n-in^dT, (2.12) 



C fl 7c fll 60 + Ar?) -T? 

where X2 — x\ + yi — y\ > 0. If we integrate with respect to ft, the integrand has a simple 
pole at jfr , which is the only pole. This pole lies outside of Cr x by Lemma 12.21 and so the 
integral with respect to £1 is zero. 

Now, suppose that the statement is true for N = K — 1. Let 

jc Rk jc Ri i 

The sum of 1(a) over <r G Bjc such that <r(if ) = K and <r ^ id is simply 

k-i 



>Cr k 



I e / c -/ c ^ nor -1 *■■■*-> 

<r(K)=K 



C K-VK-^ diK = 



by the induction hypothesis. Now we will show that 1(a) = for each a € §k such that 
a(k) = K where k 7^ K. Let Q a = {a(k + 1), • • • ,a(K)} and q = minQo-. For a fixed k, 
observe that 

i<j,k, k<j k<i<j, 

a(i)>a(j) o(i)>o(j) 



r K 



Substituting £# = ^1, so that r] € Cr where i? = Y\ { Ri in 

JCr k JCr i i<j,k, k<j k<i<j, 

< J ( i )> cr (i) <r(i)><r(j) 

n C7' w_1 ) ^- tf *- i ^-- iw " ,h " 1 ifi . ■ . #«, 



and ignoring the sign of 5-matrices, we have 

/ / ' ' ' / II ^(*Mi) II 5 <^Mi) ( 2 - 13 ) 

JC R JCr k _ 1 JC Rl i<jk ^ k<i<j, 

<r(i)>tr(j) <r(i)>ff(j) 



A^ + ieS*)-ie£s/ & W« + 



jt^ (?) 

i7^fc,er _1 (<j) 

If an S'-matrix in 

n s °(i)°(j) n 5 »(*>(j) 

i<j,k k<i<j 

°"( i )>< j 0') o"W>o"(i) 

contains £ g , the 5-matrix much be in the form of S^ (a > g) because q = mmQ a , and 

observe that the exponent of £<j CT is nonnegative. Hence, if we integrate with 

respect to £ q , there are singularities at 

v 



n 



i^K,q' 



,u+^a 



and 



where j > k, and possibly at 









from 5o,q if A ^ 0. If A = 0, there is no pole in S aq . By Lemma 12.11 Lemma 12,21 an d 
Lemma 12.31 all these poles lie outside of Cn q , and hence the integral with respect to £ g is 
zero. Therefore, 

£ w = °- 



This completes the proof. □ 

Remark 2.2. If p or q is zero, the model is totally asymmetric in the sense that particles 
move only in one direction. However, unlike the TASEP, the transition probability is not 
expressed as a determinant as long as \x ^ 1. 



3 The position of the mth leftmost particle 



Let x m (t) be the position of the mth leftmost particle at time t. In this section we find the 
probability that the mth leftmost particle is at x at time t, which is denoted by F(x m (t) = x). 



1 Notations 



Let S = {si,--- ,Sk} C N with Sj < Sj (i < j) and £ be a /c-dimensional vector in 
Define a /c-variable function on C k 



I(& si,---,s k ) = J(£; si, • • • , s k ) ■ ( ]J 6 - l) (3-14) 



where 



J& si, ■ ■ ■ , **) =- It „ + xi^- 6, ' n se5 (6 " 1) 



j <? - ■ ="-! =-j =°j .ilsegv 



for S 1 with |S| > 2. If 5 is a singleton with s £ S, then we define 

1 



J(Z;s) 



6-1' 

and so J(£; s) = 1. For the empty set, we define /(£; 0) = J(£; 0) = 1. 

Let U = {m, ■ ■ ■ ,u N } C N with m < Uj (i < j) and suppose that fu(£) = f(£ ui , ■ ■ ■ ,£ UN ), 
as a function of £ Ui , is analytic in the set {£ Ui € C : < |£ Ui | < R Ui \ for each £ Ui . For 
5 = {si,- • • ,s fe } C U let 

M£;si,--- ,sfc) = /(£ U1 ,--- ,C«jv)k Ui =i, «^s- (3.15) 

Define a map ^ : t/" ~~ ^ ^ by gu( u i) = *■ That is, <7j/ maps the zth smallest element in [/" to 
i. Define 

ses 







l 1 * J 


fjL — A 




[N]\ 


= [N][N- 


-!]• 


•[1], 


"AT 

771 


[AT]! 


[m]![iV-m]! 



for SCC/, and let U P = U \ P for P C U C N. 

Recall the notations in [TBJ for the use of the identity (1.9) in |16j . 

and 



with [0]! = 1. 



2 The position of the rightmost particle at time t 



The way of finding P(xjv(i) = x) is the same as the way of finding T{x\(t) = x) in the 
ASEP |16j . That is, we sum the transition probability over all possible configurations. Let 
xn = x and xjy-i = x — v\ — ■ ■ ■ — Vi for 1 < i < N — 1. Then the sum of Py (X; i) over all 
possible configurations is written as a multiple geometric series which is convergent 



oo oo 



(±\ N 



2-^'" 2-*< \2iriJ J Cn J Cb 

Vl =0 v N ^±=0 Cr n Lr i 

iV-1 N-k 

e *. n (n ^yais-"*- 1 ^') *■■ ■««». 

ueSjv fe=l i=l i 

Summing over each Uj, the integrand becomes 

"N-ltN-2 



e ^ — : — ;;r (2) ^-.i — n (n«.-0n(r»- 1 ^-"). 

,t£. (Wl) • • • W) ~ l)^ 1 ' ' ' £<r(N-l) -I)'" &1 - I) V V 7 V 7 



and A CT is written as 



A CT = sgn a 



Ui<j 0* + A£r(i)6x(j) - £a(j)) 



Letting £j — > ^at_j + i in the identity in (1.6) in [16] and then applying the identity to 



^2 sgnaJI^ + A^)^)-^))- 



-N-ltN-2 



V(l) V(2) •'•?<r(AT-l) 



v _ v , , (&(i)-"6t(jv) -l)(&(i) •••&(i\r-i) -!)••• (6,(1) -1) 



the sum becomes 

ll N(N-l)/2^i<j^ i ~& 
H 1 



With this result we obtain the formula for P(xjv(i) = x) of the MADM, which is to be 
compared with ¥(x\(t) = x) of the ASEP |16| . 



Proposition 3.1. With ^,Cp i in Theorem 2.4 and J(£; 1, • • • , N) given by {3. 14\ ), we have 



Mt) = x) = " / • • • / m i, • • • , n) n (£-»- v<«*) dei - <*&■ 

(3.16) 



3 The position of the leftmost particle at time t 

Lemma 3.2. Let S be a nonempty subset of U = {u\, ■ ■ ■ ,un} C N with m < Uj (i < j). 

Then, 

y^(_l)|5|+l c Sc/[5]-|5| _ ]_ 

5 

where the sum runs over all S and c is a nonzero constant. 



Proof. Let 

SJP :={S CU: \S\ = n, u x G 5} 

and 

4°) :={SCU: \S\ =n,u x i S}. 

Then, there is a bijection from Sn to 4-1 that maps S € S n to S' = S \ {u\} € 4-1- 
It is easy to check that 

(_l)l-S|+l c SdSHS| + (_l)|S"|+l c Sc/[5']-]-S"] = 

for 2 < n < N and the only remaining subset after summing over all nonempty subsets is 

S = {u\} for which 

(_l)|5|+l c Sc/[S]-|S| = L 

D 

Lemma 3.3. Let U = {u\, ■ ■ ■ ,un} C N with Ui < Uj (i < j). Given a nonempty subset 
P = {pi, ■ ■ ■ ,pk} $! U let S = {s\, ■ ■ ■ , Sk} C Up with S{ < Sj (i < j). Then, with /j, and 



10 



Cr 2 in Theorem \2. 4\ and I(£;s\,- • • , S&) and /[/(£; si, • • • ,Sk) given by flff. 1$ and Ii3.15\) , 
respectively, we have 

Y] / ••• / ( TT&i) J (^ s i'"' > s fc)/c/(C;si,--- ,s fe ,p) d£ si ■■■d£ Sk d£p 

t^o Jc Rp Jc Rsk Jc Rai vy y 

= - / ■■■/ J(£;si,--- ,Sk)fu(£;sx,--- ,s k ,p) d£ Sl ■■■d£ Sh d£ p 

Jc Rp Jc Rsk Jc Rsi 

\\ fc(fc-l)/2 



+ (2vri) A 



/i 



C fl H 



fu(C;p) d£ F 



(3.17) 



w/iere J Cj? = J Cj? • • • / Cfl and d£ p = d£ Pl • • • d£ Px . Also, f v (£; si, • • • , s fc ,p) and fu(^p) 
imply fu(£; si, • • • , s fe ,pi, • • • ,p^) and fui&Pi, • • • ,Pk), respectively. 



Proof. In the first step, we change the contour Cr s to C Ta with r Sl < (l/R Sk ) k 1 < 1 in 
the left hand side of (13,171) , Then, 

/ / "/ ( TT^si ) ^(C;si,--- >Sfc) /t/(C;si,--- ,Sfe,p) d£ Sl ■■■d£ Sk d£p 
Jc Rp Jc Rsk Jc Rsi \ =1 ' 

= ••• / (n^O J (^ s i'"' > s k) fu(£;s lr -- ,s k ,p) d£ sl ■ ■ ■ d£ Sk d£p 

JC Rp JC Rak JCr Sl i = 1 

+ (27^)-^ f f ... f 

V JC R „ JCn„ JC„„„ 



C R P JC Rs k 



Y[Csi) I(&S2,--- ,s k )f u ^;s 2 ,--- ,Sk,p) d£ S2 ■ ■ ■ d£ Sk d£, 



fc"NP 



i=2 



where the second term on the right hand side is the residue at £ ai = 1 and k _ t is from 

rr ^ s i ~ *sj 
f^.M + A^^-6 

Note that if A 7^ 0, there are other poles of £ Sl in I(£; s\, ■ ■ ■ , Sk) but these poles lie outside 
of Cr s by Lemma 12,11 Hence these poles do not produce residues. If A = 0, there are 
no poles of £ S1 in I(£; si, • • • , s^). In the second step, we change the contour Cr 3 in the 
residue obtained from the first step to C r with r S2 < (l/R Sk ) k ~ 2 < 1, which produces 
a new residue. We repeat this procedure until we obtain the integral with respect to only 
£pi>" ' >£»#-■ I n general, in the Zth step with I > 2, the contour Cr s in the residue obtained 



11 



in the I — 1 step is changed to C Ts with r Sl < (l/R Sk ) k l < 1. Thus, after the fcth step, the 
original integral is written as 

^2 c i / •••/ (n^J I ^ Sr " i s k)fu(C;sr-- ,Sk,p) d£ Sl ■■■d£ Sk d£ p 

1=1 -> C R P -> C Rs k Jc ^si i=l 

where 

'l if 1 = 1, 

ci = < 



(27ri)'- 1 






if 2<l<k. 



Now, the geometric series X)S=o ( Tli^ ?s, ) converges since | n«=/ ^ Si | < 1- In the Zth term 
where 1 < I < k, 

oo k 

^2 (II £*<) J (& ^ ■ ■ ■ > s fc) = - J (£; s r • • ,s k ), 

2=0 i=Z 

and then changing C Ts back to Cr s , the Zth term becomes 



-Q / / ••• / J{i;s r -- ,8 k )fu(€;sr-- ,Sk,p) d(, Sl ■ ■ ■ d£ Sk d£ p 

JC Rp Jc Rsk '' c Rs l 

+ c i—j—i I I ■■ I J (C; si+i ■■■ , s k )fu(b si+i ■■■ ,s k ,p) d£ Sl+1 ■ ■ ■ d£ Sk d{,p 

P JCu„ JCn.. JC R .. . . 



>C Rp JC Rsk J CRs l+] 



and the Zcth term becomes 

-Ck / J(£',8k)fu(£;sk,p) d£ Sh d£ p + c k (2m) fu(£)p) d£ p . 

Jc Rp ->c Rsk Jc Rp 

The positive term of the Zth term cancels the negative term of the (Z + l)th term for 
Z = 1, • • • , k — 1, and hence the only remaining terms are the terms in the right hand side 

of flSTrp . o 

Proposition 3.4. Let 1/2 < p < 1 and U = {u\, ■ ■ ■ ,un} C N with m < Uj (i < j). Let 

A a be given by \2.6\) with a permutation a on U. With Cr u . as in Theorem \2.J\ and fu(,0 
given by \3.15\) , we have for 2 < N € N, 



£ Q*YJc ~'J C E^lKlI^o^MOdeux-de^ 



N N 

\2niJ J Co J Cd ^ 

= E C5 (2^) / "/ J (^ s i'"' > s fc) /c/(C;si,--- ,s fe ) cZ£ sl •••<&;. 

w/iere 



5 "'" JCR s k Jc Rs 1 



*k 



\Y,u[S}-\S\ 



^u[S]-~\S\(\S\+l)/2 

and the sum runs over all nonempty subsets S = {si, • • • , s k } C U with s, < Sj (i < j) 
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Proof. We prove by induction. Let U = {u\,u 2 }. Then, 



[ [ Uuz+Suzu^fuiOd^dU 

JCn... JCn... V ' 



\1-Kl) J c J c 



I I euju{i)di ul di U2 + (^-) 2 f f s U2U1 e ul M0 dt ui d? v 

JCr„„ JCn„, X ^l) J C „ Jc r „, 



(-Y 

\2txi) jc jr„ ~ \zm/ jr Jc 

Note that the first integral in the right hand side has been obtained from changing R U2 to 
r U2 < 1 which is permissible because the integrand is analytic in the circle with radius R U2 
except at the origin, and the second integral in the right hand side has been obtained from 
changing R Ul to r Ul < 1 which is also permissible because the pole from S U2Ul lies outside 
of Cr u by Lemma 12.11 Summing over z and then changing contours back to the original 
contours, the first integral in the right hand side becomes 

(— ^ f I r-^MOd^dU (3.18) 



\2-Ki) 



Jc Ru2 Jc Rui ^-U \2ttiJ J CRui 

and the second integral becomes 

(t-)V / S U2Ul -^—fu(0 ^ Ul d^ U2 (3.19) 

JL Ru 2 JL Ru 1 SU1 r Lr u 2 

where \j \x is due to evaluating S U2Ul at £ Ul = 1. Now the sum of double integrals in the 
right hand sides of (|3.18p and (|3.19p is 

\(__\ I / Wjl ~SM2J(sm?M2 ~ 1) f (C\ AC AC 

X W0 J Cru2 J Crui fr + xuu - U)(U - i){U - i) M5j Ul u ' 

Hence we have 

+ (^~) I ^(e;«i)MS;ui)de« 1 + (-)(2^) / i(^u 2 )fu(^u 2 ) du 

where /(£; u\) = I(£; u 2 ) = 1. 

Suppose that the statement is true for N = K — 1. Pick up p € {ui, • ■ ■ , un}- We first 
sum over all a with cr(l) = p, that is, we compute 

oo N N 

£ (^r /, -I E ^ II (n^^^/^e) ^ • • • «%,„. (3.20) 

2l,-,2jV-l=0 R ™N JC Ru 1 a with fc=2 i=fc 

<r(l)=P 
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Let us write A„ as 



ff(j)<p, 2<i<j, 

i>2 <t(*)>»0') 



Changing the order of integration, which is permissible by Fubini theorem, we write (J3.2 1 



as 



1 r °° 

2^?/.L ^ 



c % zi=0 

N N 



^ \ 9-7T7 

Z2,— ,*JV_1=0 



2vri/ ./ c y Cn 



e n * 



er{i)ff(J) 



(2) <t with 2<?<j, 
°"( 1 )=Pcr(i)>o-(j) 



n(n^) 

fc=3 i=fc 



^fc-i 



n Va) ( n ^o fr(£) d ^(2) ■ ■ ■ ^w 



where 



a(j)<p 
J>2 



n vo)(n^) m& 

'"(j)<p. u i+v 

i>2 



*; P , 



as a function of £ CT (2), ■ ■ ■ ,£,cr(N)> is analytic in the set {£<r(i) € C : < |£ ff ft)| < -R CT (j)} for 
each i = 2, • • • , N by Lemma I2.ll and the assumption of fu(0- Hence, by the induction 
hypothesis, we have 



^^flp 2l =0 5c(7p -"-^fc l/c «»i - " 



a<p, 



a<p, 



i£U P \S 



(3.21) 



where S = {s\, • • ■ , s^} with Sj < Sj (i < j) and P = {p}. Note that eg in (|3.2ip is defined 

for S C Up. If we denote the number of elements in the set {a : a < p and a £ S} by I, 

then 

r -i-r 1 /AnsuCp)- 1 -' 

L 11 Sn 



'pot 



a<p, 



£a=l ^ 



Hence (|3.2ip is written as 



E/ \\9u(p)-i--l / 1 \fc+l ^ 
c 4- u- L 



5cC/p 



/" 



2m 



Cr*, J CRs,„ Jc r 



Z1 =0"^R P "^Rs k 



I(€)8l,--- )Sfc) JJ •S'pa 



a<p, 



i&Jp\S 



(3.22) 



«gs 
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Now, we can apply Lemma [3T3l to (|3.22p since Y\ a <p aes Spa[fu(Q] &=l, satisfies the 

ieu P \s 

condition for the /[/ in Lemma I3T31 Observing that Y\ a< P> Spa = (VaOS we have, 



aGS 



e«=i 



after applying Lemma 1331 

]J <%** A/(£)l &=i, d 6i • • • d£s k d£ p 

a<p, aeS i£U P \S 

In the second sum, 

/Ayt/(p)-i/l\fc(fc-i)/2 _ | 5 | +1 /A\=Ui»['S(l-l'S f l+si/(p)-i 

and by Lemma |3, 21 the second sum becomes 

. >5l/(p)-l-2 
Now we express eg I ^ ) in the first sum in terms of 5 := P U S C U. 



Ma) = ( - 1} " 



(-l) fc — r 



•1 



-A' 
V- 1 

\T$ 9u P (*i)+(k-l)+9u(p)-2k-l 

,fc _^ 

Ei^p( s i)+( fe -0+ffa(p)-fe-fe(fe+i)/2-i 



\Ec/[S]-2|S|+l 

= f-iyg|-i A -=c- 

V ; ^ v [S]-\S\(\S\+l)/2 ■ S- 

Hence, ()3.23j) is written as a single sum 



, SfcJ 



scc/p "" J CR v Jc Rsk JC Rsi 



II 3«M£)I &=i, <&i ■■•<&*#* (3-25) 

a<p, aGS ieU P \S 
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where S is any subset of Up. When S is empty, JIcKp aes Spa is defined to be 1 and the 
integral is over only £ p . Now, we sum f)3.25j) over p G U. Observe that for a nonempty 
S CU P with \S\ > 2, 

jft. e a \ TT <? - U.i<j(€sj ~6j) . ,jrT (m_+ Ag s Cp-6) 

( j ii, V " IL^O* + A6.6, - 6,) IL eS (& - 1) ( } /J, (m + A6C P - Q 

aes ses 

= J(^;ai,--- ,a fc ,p)(-l) ^ 77 TTYpj 77 — (£ p - 1) 

S<p S>p 



Uses (£p ~~ £s 



J& si,-'-, Sk,p V ise ^ \,Z, SJ (£ P - I)- (3.26) 



Note that if \S\ = 1, the last equality in (|3,26|) is still obtained. Hence, the sum of (|3.25p 
over p € U is written as 

s cU \s\=k-i 



C§ \2m) „ 

n 

sGS, P eS\S 


f ...J j 

'c Rsk Jc Rsi 

ix + xtsip - is 


sp ~~ ss 



ieu\s 

where S = {si, • • • , s^} with fe > 2 and S C S C U. Using the identity (1.9) in [16] for the 
sum over all S with \S\ = k — 1 in S, we finally obtain 



E^l^i) / "/ J (^ s l'"' > s k)fu(£;si,--- ,s k ) d£ Sl ■■■d£s 

sen "" Rs k Rs l 



where 

c, := (-1)I 5 I +1 . , _ , 

If S 1 is empty, c~a = ca in (|3.25j) . and (|3.25[) becomes 



" S \2m/ j c 



c 

1 \ ''71-// /. 



)y fu&p)dt P 



for each p. This is combined with (|3.27p to complete the proof. □ 

Corollary 3.5. Let 1/2 < fi < 1 and S 1 = {si, • • • , s^} be a nonempty subset of {1, • • • , N} 
with Si < Sj (i < j) and N > 2. In the MADM with N particles the probability that the 
leftmost particle is at x at time t is 



F( Xl (t) = x) = Y,cs(^-) k I ■■ j J(e; *!,•••, *fc) II (^~ w -~ 1 e c(e ' ) *) deax • • • «-* 

(3.28) 



ic Rsk Jc Rsi sgS 
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where 

CO = <-l)\ S \ +1 A 

and Ti[S] is the sum of all elements in S. 



Proof. Choose U = {1, ■ • • ,N} and f v {0 = Ui^i W V^* in Proposition E3J Then 
P(xi(i) = x) is immediately obtained. □ 



4 Main results : the position of the mth leftmost particle at time t 



Theorem 3.6. Let 1/2 < /x < 1 and U = {u\,--- ,un} C N with m < Uj (i < j) for 

N > 3. Let A„ be given by Ii2.6\) with a permutation a on U . With Cr u , as in Theorem \2. 
and /[/(£) given by k°3.15\) . we have for 1 < m < N 



m—l m~k 

e e (£)"/. ■■•/. E^n(n«*«> 



^ \2tHJ J Cn Jc„ z 

Zl,— ,z N - m =0vi,--- ,tl m _l=0 """AT H «l cr&N k=l i=l 



N N 

Zk- 



k=m+l i=k 



^ C5 V2vri) / '"/ J (£'' Sl '"" > s k)fu(&sx,--- ,s k ) d£ sl ---d£ Sk 



\S\>m JCR s h Cr *i 



where 



cs = (-l) |A|+m (/iA) 



\S\+m (llX \m(m-l)/2 



\S\-1 
\S\ — m 



X T,u[S]-m\S\ 



^u[S]-\S\(\S\+l)/2 

and the sum runs over all S = {s±, • • • , Sk} C U with Si < Sj (i < j) and k > m. 



Proof. Let Q C U such that \Q\ = N — m, and P = U \Q = {p±, ■ ■ ■ ,p m }- Let q be a 
permutation on Q, and let us denote by S q a set of a with <r(m + i) = q(i) so that 



for a € S q . Observe that 



Q = {o-(m+ !),••• ,<t(JV)} 



E = E E E 
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where the first sum in the right hand side is over all partitions {P, Q} of U with \P\ = m 
and \Q\ = N — m. We first sum over all permutations in S q . Observe that geometric series 
over Vi converge for all i and A a is written as 

k= n s ff ({)<r(j) n s a{i)u{j) w s a{i)u{j) . 

i<j<m, i<ni'<j> m<i<j, 

<r(i)><7(j) o"(i)>cr(i) a(i)>a(j) 

Since a(m + i) is fixed to q(i) for each a € S q , we have, after summing over all Vi, 

£ (s)"j£ -L 

V(l) V(2) ' ' 'Sa(m-l) 



y TT s . . Ml) M2) 

aeS q i<j<m, ata3 (&M ' ' ' &(m-l) ~ X ) (&(1) ' ' ' 6r(ro-2) - l) • • • (&r(l) ~ l) 

°"( i )> CT (i) 

AT AT 

II S a(i)a(j) [[ S a(i)(T(j) 11 (II^WJ fu(0 d£ Ul ---d£ UN , 

i<m<j, m<i<j, k=m+l i=k 

<r(i)>o{j) °"(*)>°"(i) 

where the sum over a € <S g is equal to 

/x m(m - 1)/2 IK;c7(l),---,a(m)) 

by the identity (1.6) in [16]. Now, we sum over q E S^-m- Let us write the sum of (|3.29p 
over q € §.v_ m as 

^-iw. E E (_L) 

Zl=0z 2 ,",2jV- m =0 Cfl «iV Cfl «l 

AT AT 

(?GSjv_ m m<i<j, k=m+2 i=k 

°"( i )>°"(i) 

AT— m 

^(^Pl)--- ,Pm) 11 \Ca(m+i)) fu(Q [[ S(r(i)<r(j) ^«i " ' ^ w w • 

j=l i<m<j, 

Observe that 

7V-m 
^(C;Pl,"- ,Pm) J| (? CT (m+i)) M£) II S v(i)<T(j) 

°"( i )> CT (i) 

as a function of ^( m+ i),--- ^(jv), is analytic in the set {^(m+i) G C : < |£r( TO +i)l < 
-R CT ( m+ j)} for each ^o-( m+ j) (i = 1, • • • , N — m) by the assumption on /[/(C) an d Lemma [27X1 
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Hence, we may apply Proposition 13.41 to the multiple integral over £ CT ( m+ i), ■ ■ ■ , £<j-(jv)- That 
is, if we change the order of integration which is permissible by Fubini's theorem, we have 



/<' 



i(m-l)/2( \ \ N 



\27TiJ J c J c 



E 



Pm n P\ 2i=0 

N N 



t L -L 



Zl 



N—m 

Yj II S <r(i)<r(j) II (n^(0) " mj (^Pl>"" >Pm) II (^r(m+*)J " MO 

(jeS]V-m m<i<j, k=m+2 i=k i=X 



] J S<r(i)<T<j) ^ct(AT) • • • d£o-(m+l) 



i<m<j, 
°"(*)>°"0') 



"Spi ■ ■ ■ "SPm I 



and applying Proposition 13.41 we have 

-i(m-l)/2 ( J_ \ m+fe 



A* 



V27ri; y Cfl 



Cu P1 ^ 1=05cQ ^C flafe JC flsi 



Or, 



i<m<j, i<m<j, 

a(i)>ff(j)eS <r(i)>a(j)eQ\S 

i(bpi, ■■■ ,Pm) [fu(0] e<=i. d &i ■ ■ ■ d&fcd&i • • • ^P« 

i£Q\5 



&Ki) : 



..ne- 



sGS 



(3.30) 



where 5 = {si, • • • , s^} with Sj < Sj (i < j). Note that eg in (|3.30p is defined for S C Q. 
For p £ P, let us denote the number of elements in the set {a : a < p, a € P} by /^ and 
the number of elements in the set {a : a < p, a £ S C Q} by l~j. The number of elements 
in the set {a : a < p, a £ Q} = gu{p) — lp and hence 



n ^(^(j) 

i<m<j, 
*(i)>a(j)£Q\S 



t*U) =1 \fi 






■iZ-iS) 



Hence, (|3.30p becomes 



/' 



i(m-i)/2(J_\ m+k sr^ c /A\E P€ i'(^(p)- J p- J p) ^ 



K27TiJ ^Q "''" 



Zl =0 JC R Pm JC Rpi JCR s k JC R S1 

I(£; sx, ■ ■ ■ , Sk) Yl S °(iMJ) II & ffaPi' " ' '^ m ) [MO] &=i, <%si ■ ■ ■ d £s k d£ Pl ■ ■ ■ d£ Pri 

i<m<j, s£S ieQ\S 

a{i)>a(j)eS 

Now, we may apply Lemma 13.31 since 

i<m<j, i£Q\S 

<r(i)><r(j)€S 
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satisfies the condition for the fu in Lemma 13.31 Hence, observing that 



i<m<j, 
<r(i)>ff(/)eS 



^(j) : 



.. - 



A^ E P gp 'p 



we have 



/' 



i(m-i)/2(J_\ m+k sp /A\E p£ p (gv(p)-i?-i§) 



VW sTq v "' 



JC Rn™ JC Rvi JC R S ,. JC Rs, 



91,"' i S k) 



J] S <r{i)<r(j) I{£;P1, '" ,Pm) [fu(0] &=1, Csi ■ ■ ■ d£ Sk d£ pi ■ ■ ■ d£ p 

i<m<j, i&Q\S 

cr(i)>ff(j)eS 



+ /i' 



\2TriJ '— -* v/i 






A \ E p£ p (ft/(p)-^) / 1 \ fc(*-i)/2 



/' 



/ ■■■ / J(6pi,--- ,Pm)[^(0]&=i, dZpi'-dtpm 

JC Rvm JC Rn iGQ 



(3.31) 



by Lemma 13.31 In the second sum of (|3.31|) . 



- A \ E pG p {guto-i?) / 1 \ fc(*-i)/2 



1 /A\Eg[SHS|+E peJ ,(sir(p)-J!') 



// 



and the second sum is written as 



»"iM a O 



1 \ m /A\E pS p (au{p)-ip) 



c Rpm j Cr Pi 



i(€;pi,'~ ,Pm)[fu(€)]ti=i, d£pi'-d£pr, 



, T ") V / A \E, e p(flu(p)-J|'-^) 

by Lemma 13.21 JNow we express — cs\ 



in the first sum of (|3.3ip in 
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terms of S := P U S. Observing that J2 p£P lp = m(m + l)/2, we have 

a S q [5H5|+E p£ p (str(p)-J?-j£) 



'A\E p6P (9!/(p)-' P p -i P s ) 



»ISI. 



E Q [S]-|5|(|S|+l)/2+E pe p(sc/(p)-^-i#) 



/' 



\|S|. 



/' 



.1)151. 



/' 



A E- ffo(si)-fc+E peP (sp(p)-£-i?) 

E- 9Q(^)-fe(fc+l)/2+E p£ P ( S p(p)-^-^) 

A E* 9Q(^)+(™A:-E pe p ^)-("i+i)fc+E pe p (su(p)-£) 
E*9Q(si)+(mfe-E pe pif)-mfc-fc(fe+l)/2+E pe p(9u(p)-i|') 



.1)|5| 



\Ei 9q(^)+E pS p ffU-(p)+Ep e p(fe-i|)-("i+l)(fe+"i)+"i(»"+l)/2 



_-]\[S|-m 



^E* SQ(*i)+E pe p 9y(p)+E pe p(fc-^)-(™+fe)(™+fc+i)/2 

A Su[S]-(m+l)|S|+m(m+l)/2 



^[sj-ISKIsi+i)^ 



Hence, (|3.3ip is written as a single sum 

Mm - 1)/2 (^-f l E"s[ ■••/ [ ■■■[ J(t;*i,- 



/i 



ScQ 



>c Rpi Jc Rsk 



I] S °(i)<r(J) fibPl' " ' iPm) [fu(0] 6=1, d ^si ' ' ' ^Sfc^; 

i<m<j, i€Q\S 

a{i)>a{j)&S 



VI 



<% Pm (3.32) 



where S is any subset of Q. When S is empty, Y\ i<m<j, S<r(i)<r(j) * s defined to be 1 and 

<r(i) ><r(j)'g S 

the integral is only over £ pi , ■ ■ ■ , £ Pm . Now, we sum (|3.32p over all partitions {P, Q} of U. 
Notice that if S is nonempty, then 



n 



s, 



<r(i)ff(j) 



1 I P s ' 



i<m<j, 
a(i)>a(j)£S 



s<p, 

p&p,sgs 
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and 



J{£;s!,--- ,s k ) I(£;pi,--- ,p m ) Y\_ S ps 

s<p, 

p&p,sgS 

n,<j(6i -6,-) n Pi <Pj (^ ~ ^i ) ( ripsp 6 ~ x ) 



(_l)E pe p'p T~T 



s<p, 

pgp,s£5 



(m + Hs€p - Zp) 



= J(C;si,--- ,s fc ,pi,--- ,p m )(-i) E ^ p/p (n^" 1 

peP 

n ^<p, (M+A6^-6)n *>P: (m+a6c p -6) 

peP,st=s peP,seS 

n s<p, (6 - Cp) n *>P: (c P - 6) 

peP,s&s peP,s£S 

= J{^s 1} --- ,s k ,pi,--- ,Pm) ff Tc T\ ( 1 Hp~ 1 i- 

llseS, P ep(£p - Zs) \£* / 

Hence the sum of (|3.32p over all partitions {P, Q} of U is written as 

m(m- 



/' 



-""£ JE ^)7 -1 ^;.,-.«) 



n 






n Cp-i)MOi&=i,.de-i---de. 



where 5 = {si, • • • , s^} with fc > m + 1 and S C S C U. By using the identity (1.9) in [16 
for the sum over all S with IS 1 ) = k — m, we obtain 



^2 c s{^—) / ■••/ J (^ s l> - -- ,Sk)fu(€; «!,••• ,Sfc) ^«i"'^«* 



(3.33) 



scu 



>c Rsh jc Rsi 



where 



(—l\\S\+m^y.m(m-i)/2 



'fc-1' 

k — m 



X T,u[S]-m\S\ 

^[si-iskisi+i)^" 



If S is empty, then /i"^™- 1 )/ 2 ^ = c § in (^321) . and ([332]) becomes 

c s(^) / '"/ j (£;pi>--- ,Pm) fu(£;pi,--- ,p m ) d£ P i---d§ 

for each P. This is combined with ()3.33p to complete the proof. 



□ 
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Theorem 3.7. Let 1/2 < fi < 1 and S = {si, • • • , s^} be a nonempty subset of {1, • • • , N} 
with S{ < Sj (i < j) and N > 2. In the MADM with N particles the probability that the 
mth leftmost particle is at x at time t is 



l (t) = X )= Y,c s (^-) k f ■■■[ msi,--- ^^nes-^-'e 6 ^ du---du 

(3.34) 



\S\>m " J ^ Rs k ° R n ses 



where 



Cs = (_l)|5|+m^ A )m(m-l)/2 

and Ti[S] is the sum of all elements in S. 



\S\-1 
\S\ — m 



X nS]-m\S\ 
^{S]-\S\(\S\+l)/2 



Proof. When m = 1 and m = N, (|3.34|) is equal to (|3.16|) and (|3.28|) . respectively. When 
3 < N, choose U = {1, • • • ,N} and fu(£) = fj. ^"^"^^fe)* in Theorem El for 1 < m < 

N. D 



Let us denote the number of particles in (— oo,x] at time t by T(x,t), which is called the 
total current |20j . The position of the mth leftmost particle at time t is related to the total 
current. That is, from the fact that 

{T(x,t) = m} = {x m (t) < x, x m+ \(t) > x} 

it can be shown f20j that 

P(T(x, t)<m-l) = l- P(a? m (t) < x). 

Hence summing (|3.34p over x from — oo to x to obtain ¥(x m (t) < x), we have 

F(T(x,t) <m-l) (3.35) 

= i-E ^^-)"/ •••/ J&su---,s k )i[es- ys e £ ^du---du- 

\S\>m m Cfls fc Cr °i seS 



4 The two-sided PushASEP 



The MADM is closely related to the two-sided PushASEP. In the one-sided PushASEP the 
pushing effect is totally asymmetric and it is possible to express the transition probability as 
a determinant of contour integrals [4] as in the TASEP. In the two-sided PushASEP, unlike 
the one-sided PushASEP, the pushing effect is partially asymmetric, that is, the pushing 
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effect is allowed in both directions. The dynamics of the two-sided PushASEP with N 
particles is as follows [2]: each particle is equipped with 2N Poisson clocks with rates pr n 
and ql n (p + q = 1,1 < n < N), where r n and l n are given by (|1.1|) and (|1.2p . respectively. 
All clocks are independent and each site can be occupied by at most one particle. If the 
nearest empty site on the right (left) is x + n (x — n) when a clock with rate pr n (ql n ) of 
a particle at x rings, then the particle at x jumps to x + n [x — n). Otherwise, nothing 
happens and the clock resumes. As mentioned in [3] the jump to the right (left) of a particle 
to the nearest right (left) empty site is also interpreted as the particle's pushing all its right 
(left) neighboring particles by one if they prevent the particle from jumping to the right 
(left). The transition probability of the two-sided PushASEP was treated in [7], however, 
its proof is unfortunately false because the proof used Lemma 2.2 in [16] which turned out 
to be falsqj. In this section we obtain the transition probability and the probability that 
the mth leftmost particle is at x at time t for the two-sided PushASEP as corollaries to 
Theorem 12.41 and Theorem 13.71 respectively. In the PushASEP the physical region is 

{(xi, • • • , xn) 6^ : x\ < ■ ■ ■ < xn} 

because of the exclusion property. The S- matrix of the two-sided PushASEP [2 , 7J is defined 
to be 

b ^- &' » + >£& -Zp ^ $a 

where Sp a is the S'-matrix of the MADM and the Bethe ansatz solution is 

N 

z^ ^llQi) 



with 



Recalling Lemma 3 in [7] 



M _ TT ^W . a 

°"( i )> CT (i) 

°"(0>< J 0') 
we have the following transition probability of the two-sided PushASEP. 

Corollary 4.1. Let A + jx = 1, 2 < /" — 1 an< ^ ^Ri (* = l; --- >N) be a circle oriented 
counterclockwise, centered at with radius Ri. Assume that 1 < R\ < • • • < Rn < c where 

I 00 if A = 0. 



The authors corrected the error in the erratum [16] 
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The transition probability of the two-sided PushASEP with N > 2 particles is 

N 

«-<**>-£ (as) / - ^n^S"""" 1 ^'') *■■■«* (4.36) 



crgSjv 



/C fljv JC fll 



^nfc7 CT(i) ~ 1+<TW " le£(6)i 



Proof. The integrand of (|4.36p is equal to 

AT 

T(i) 

by Lemma 3 in [7] where vlo- is the coefficient of the Bethe ansatz solution of the MADM. 
Since the only difference between (|2.9p and (|4.36p is the exponent of ^ c u\ , we simply modify 
the proof for the MADM. The exponent of £i in (|2.12p should be replaced by X2 — X\ + 
2/2 — Vi — 2 which is nonnegative in the physical region of the two-sided PushASEP, and the 
exponent of £ g in ()2.13[) should be replaced by 

x ff -i(g) - x k + y K - y g + q - o-' 1 ^) - K + k 

which is also nonnegative since x ff -ir q \ — x^ > a^ 1 (q) — k and yx — y q > K — q in the 
physical region of the two-sided PushASEP. The rest of the proof is followed by the same 
argument as in the proof of the MADM. □ 



From the alternate form of (|4.36p 

N 



(s*)7 e -j e ^n(c„r (, "" , "' < ' ,) " 1 ^ w ' 



*<*<>- x; =j . ■■• . ^nc ( ,r"*^ «.•••«»• 



it is confirmed that the configuration (xj)j in the two-sided PushASEP is mapped to the 
configuration (xj — i)j in the MADM. Hence, we immediately obtain ¥(x m (t) = x) for the 
two-sided PushASEP by this mapping and (|3.34p . That is, in the two-sided PushASEP, 

P(x m (t) = x) (4.37) 

= E «(r)' / ■■■/ J tf; «i.-- ■.«*)IIff" (I ""' ) " v(W ' *i •••*•*• 

life vw 7c ^ fe 7c ^ ii 

In the totally asymmetric limit A — Y 0, only the subset S = {1, • • ■ , m} that satisfies 

m(m-l) 

h X[6J — m|6| = 
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survives in (]4.37p. Hence, in this limit, we have 






F(x m (t) = x) 






= (i)7 .../ , (ftl ,.. 

™™J Jc Rm Jc Rl 


m)nc (2/i_i)_i e £(cs) *^i- 


• • d£ m . 



This result tells us that the probability for the mth leftmost particle's position is deter- 
mined by initial positions of the only first m particles. This makes sense physically because 
in the limit A — > 0, particles are governed by the TASEP dynamics in the left direction but 
freely jump to the right neighboring site by pushing other particles if they occupy right sites. 
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